Ergodic properties of a simple deterministic 
traffic flow model re(al)visited 



Michael Blank** 
June, 2002 



Abstract. We study statistical properties of a family of maps acting in the space of inte- 
ger valued sequences, which model dynamics of simple deterministic traffic flows. We obtain 
asymptotic (as time goes to infinity) properties of trajectories of those maps corresponding to 
arbitrary initial configurations in terms of statistics of densities of various patterns and describe 
weak attractors of these systems and the rate of convergence to them. Previously only the 
so called regular initial configurations (having a density with only finite fluctuations of partial 
sums around it) in the case of a slow particles model (with the maximal velocity 1) have been 
studied rigorously. Applying ideas borrowed from substitution dynamics we are able to reduce 
the analysis of the traffic flow models corresponding to the multi-lane traffic and to the flow with 
fast particles (with velocities greater than 1) to the simplest case of the flow with the one-lane 
traffic and slow particles, where the crucial technical step is the derivation of the exact life-time 
for a given cluster of particles. Applications to the optimal redirection of the multi-lane traffic 
flow are discussed as well. 

Keywords: dynamical system, traffic flow, substitution dynamics, attractor, rate of conver- 
gence, large deviations. 
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1 Introduction 

Let Xm :={# = (... x^\XqX\ . . .) : X{ £ Am, i £ 2Z} be the space of be-infinite sequences 
(which we also call configurations) from the alphabet Am = {0,1,2,... , M}. We equip this 
space with the metric 

oo 

dist M (x,y):= ^ ( M + l) -|<l k» - Vi\ 

i=—oo 

and consider a map Ti m '■ Xm — * Xm from this metric space into itself: 

(T ltM x)i := Xi + min{xi_i, M - Xi} - minjxj, M - x i+1 }. (1.1) 

One can interpret the i-th coordinate of x £ Xm as Xi particles and M — Xi holes (empty places) 
located at the site i of the integer lattice 2Z. Then this map can be considered as a discrete time 
/ discrete space model for multi-lane highway traffic when a particle (vehicle) at site i of the 
lane j can switch to any other lane j' (nonnecessary neighboring) whenever it does not disrupt 
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the motion of other particles, i.e. the sites i, i + 1 of the lane j' are not occupied. From the 
point of view of probability theory the dynamics of this map is a deterministic version of an 
asymmetric exclusion process, i.e. the motion of a collection of random walkers constrained to the 
nonintersection assumption (see e.g. |], |6|). Traffic flow phenomena have attracted considerable 
interest during last decade both from the applied and theoretical points of view. For the general 
account on these matters we refer the reader to recent reviews |4| (and numerous references 
cited there) and in this paper we shall concentrate only on the mathematical background of 
deterministic models of traffic flows. 

We shall refer to the system (7i m,-Xm) as the slow particles model, and to take into account 
traffic flows where particles can move with the (maximal) velocity v > 1 (a fast particles model) 
we consider a family of maps T Vj m '■ Xm - > Xm describing the M-lane traffic flow model with 
the maximal velocity \v\, i.e. a particle in this flow can move to the right (left if v < 0) by at 
most |u| positions if those positions are not occupied. 

To simplify the notation we shall drop the indices if they are equal to 1, i.e. T<i means the case 
v = 2, M = 1 and T\ 3 means the case v = 1, M = 3, while T stands for the case v = M = 1. 
By a dual configuration for the configuration x G Xm we mean a configuration x* G Xm such 
that x* = M — Xi Vi. The operation of taking a dual can be applied also to the map by means 
of the relation: T* M x* := (T V; m%)* Vx G Xm- 

To illustrate the usage of the dual operation consider the slow particles model with 'smart 
drivers', who anticipating the motion of at most m cars ahead, may move to an occupied site 
ahead of it with the maximal velocity 1. Example for the case m = 2: (OHIO) —* (01011), where 
(•) denotes the main period of a (space) periodic configuration. It is straightforward to show 
that this model is described by the map Tx = (T^ m< ix*)* = T* m x x. 

By a word A we shall call any (finite or infinite) sequence of elements ai G Am and introduce 
the notion of the density of a finite word A in a finite word B as 



p(B,A) 



1 

W\ 



\B\-\A\ 



mm' 



i=l 



B 



i-l+j 

A, 



;i.2) 



where \A\ is the length of the word A, Aj G Am is the j-th element of the word A, [a\ is the 
integer part of the number a, and we set 0/0 = 1 here and in the sequel. In the case M = 1 the 
number p(A, B) G [0, 1] and is equal to the number of occurencies of the subword B in A divided 
by the length of A, while in the general case M > 1 we have p(A, B) G [0, M] and the formula 
( |1.2j ) takes into account multiplicities of those occurencies. Example: p(255,12) = £(2+2) = 4/3. 

The generalization of this notion for an infinite word / configuration x G Xm leads to the notion 
of lower/upper density: 



p±(x,A) :-- 



lim 

n,m—> 00 



sup 
inf 



p(x[-n,m],A), 



where (and in the sequel) lim sup corresponds to the index + and lim inf to the index — , and 
x[n, m] a subword of the word x which starts from the position n and goes till the position m in 
the original word. The asymmetry with respect to n and m is necessary to take into account the 
possibility to have left and right 'tails' with different statistics: for x := . . . 00001111 ... we have 
p~(x, 1) = and p + (x,l) = 1> while p(x[—n,n],B) n _z!°? 1/2. Observe also that for a (space) 
periodic configuration (A) := . . . AAA ... we have p-({A), B) = p + ((A),B) = p((A),B) for any 
pair of finite words A, B. 

For a system of particles on a lattice one can define its average velocity as follows. First, for 
each particle in a configuration x G X we define its 'local' velocity as a distance by which it will 
move on the next step of the dynamics: 



V(x,i) := mm{v, x^, min{j — i : j > i,Xj = 1}}, 
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Figure 1: Fundamental diagram for T Vj m- dependence of the average velocity V or the flux 
on the density p = p(-, 1). 

and, since for M > 1 a site i in the configuration a; € 1m may contain several particles (i.e. 
Xi > 1), we sum up their velocities to get V(x,i). For example, in the case of the map T^m 
we have V(x,i) := min{xj,M — Xj+i}. Note that the 'local' velocity is well defined for each 
site i of a configuration x G 1m (independently on the presense of a particle there), indeed, if 
Xi = (i.e. there is no particle at this site) we have V(x, i) = 0. Now we define the lower/upper 
average velocity as 

V±(x) := lim [ S T ] — r -i- V V(x,i). 

n,m-^oo ^ mf J p(x[-n,m], 1) • (n + m + 1) 

Often it is more suitable to work with another statistics, called flux, equal to the number of parti- 
cles crossing a given position on the lattice per unit time, i.e. <3>(x[— n, m]) := n+ l n+l Y^L-n ^( x > *) 
Thus we define the upper/lower average flux as 



$+(x) ■= lim ( SUp ) - V V(x i) 

n,rn^oo ^ inf / n + 7X1 + 1 



n + m + 

t— — 11 

We shall use also the notation to indicate the maximum velocity if needed, and 0j := 00 . . . . 

i 

The connection of the flux to the densities is given by the following simple result. 

Lemma 1.1 <&^{x) = Ya=i P±{ x -> 10* ) f or x £ X, in particular ^± (x) := p±(x, 10). 

Proof. By definition we have 

*%>(x) = v ■ p±{x, 10„) + (v - 1) ■ (p±(x, 10„_i) - p±(x, 10„)) + . . . + 1 • p±(x, 10) 

= («-(«- 1)) • p±(x, 10„) + ((« - 1) - (« - 2))p±(x, 10„_i) + . . . + 1 • p±(x, 10) 

v 

= ^2p±{x, 10j). 
i=i 

□ 

The main results of the paper are the following statements. 

Theorem 1.1 (Invariance of densities) p±(T^ M x, A) = p±(x,A) for all x G 1m «^ t G 2Z + 
i/ and only if A G {0, 1}. 

Denote by Free„ := {x G 1m : V(a;,i) = v ■ x% Vi G 2Z} the subset of configurations 
where all particles have the maximal available velocity and thus move independently. Clearly, 
T„.M(Free„) = Freet, and T„.M(Free*) = Free*. It is of interest that Free n Free* / 0. 
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Theorem 1.2 (Convergence) The set Free„UFree* is the only locally maximal weak attractor of 

the dynamical system (T v ^m,Xm), and for x G Xm we have T* M x i " 'f ^ + ' | ~ "j^ 1 

Irree^ ij p—{x, L) .> 



u+1 



Theorem 1.3 (Limit flux) ^ Km x ) ^ (p±(*, 1)), tofcere F„, M (0 := _ ^ fj^^ 1 ■ 

Denote by p p a product (Bernoulli) measure with the density pM on the space of sequences 
Xm- 

Theorem 1.4 (Typical dynamics) For p p -a.a. x G Xm we have p(x, 1) = pM and distjv/ (Xj j^x, PreeU 
Free*) < Af"*/^ 1 and limsup^ £ £™ = _ n V^m*, 1) = F v>M {p(x, 1)) /or any 7 G (0, 1)'. 

In this result one can use instead of p p any probabilistic translation invariant measure with 
fast enough decay of correlations (see Lemma 2.11| ). 



Proofs of Theorems |l.l| -p~4| are based on the reduction of the general case v, M > 1 to the 
simplest one v = M = 1. For v = 1, M > 1 this reduction boils down to the proof that a 



multi lane traffic flow can be represented by a direct product of one-lane flows (see Theorem 4.1 
describing the 'sawtooth redirection' construction). In the case v > 1,M = 1 we make use of 
a specially constructed substitution dynamics (see Lemma |3.1| ) to prove the reduction, while in 
the general case v,M > 1 we combine these two arguments. The main technical step of the 
analysis in the case v = M = 1 is the derivation of the exact life-time for a given cluster of 



particles, i.e. the number of iterations after which it will disappear, described in Lemma 2.4. 
Note that earlier only very weak (and unnaturally large) estimates of the life-time type were 
known (see, e.g. [@, @, §). 

We provide also the analysis of the rate of convergence to the limit of various statistics for 
space periodic, regular, and 'typical' initial configurations based on large deviations estimates 
(see Section |||) , and study the dynamics of a passive tracer in the flow of fast particles (Section 0) . 

It is clear that the main problem in the study of traffic flows is the analysis of 'traffic jams' 
(without them the dynamics is trivial): we shall say that a segment x[n, m] with m > n cor- 
responds to the jammed cluster if max{V (x,i)/xi : n < i < m} < v and V(x,n — v)/x n — v = 
V(x,m)/x m = v, i.e. all particles inside of this segment do not have the maximum available 
velocity. Note that in the case v = M = 1 the jammed cluster is the same as the cluster of 
particles. 



2 Dynamics of slow particles (T, X) 

This model has been introduced originally in [|ll|] for the case of a traffic flow on a finite lattice 
(say of size L) with periodic boundary conditions and studied numerically in a large number of 
publications. It is straightforward to show that this case corresponds to the restriction of the 
map T to (space) L-periodic configurations. The first 'quasi'-analytic result for the L-periodic 
case has been obtained in M for 'typical' initial configurations of length L. However the first 
complete proof appeared only in [0], where regular initial configurations on the infinite lattice 
were considered as well. In this section we shall study the problem for all initial configurations, 
using a rather different and more simple approach than the one in J|]. 

Let us start from the analysis of lower and upper densities. Note that if the lower density 
coincides with the upper one, i.e. the limit value exists, we cal this common value the density 
p(v)- Example when they do not coincide: 1 = p + (. .. 111000 1) ^ .. 111000 1) = 0. 

Lemma 2.1 p±(x, 1) = 1 — p^(x, 0), and thus p±(x*, 1) = 1 — p^(x, 1). 
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Proof. By the definition of the lower density we have 



P-(x, 1) = liminfp(a;[— n, m], 1) = 1 — limsupp(a;[— n, m], 0) = 1 — p+(x, 0), 

since p(A, 1) • |A| + /?(^4, 0) • |A| = |A| for any finite binary word A. The derivation for the 
upper density follows the same argument, while the second statement follows from the identity: 
p(x*[— n, m], 1) = p(x[— n, m], 0) = 1 — p(x[— n, m], 1). □ 



Lemma 2.2 ^(x,^) > p±(x,B) ■ p(B,A) for any configuration x E X and any pair of finite 
words A, B. 



Proof. If A^-B the inequality becomes trivial, since p(B, A) = 0, while p±(x, A) > 0. Assume 
now that A C B. Then Vn, m £ TL we have p(x[— rt, m], A) > p(x[— n, m], B) ■ p(B,A) because 
the right hand side takes into account only those enclosures of B to x when the word B belongs 
to a segment z[i,j] = A, while there might be other enclosures as well. □ 



Proof of Theorem 1.1 in the case v = M = 1. Let us prove first that p±(x, 1) = p±(Tx, 1) for 
all x € X. For any n, m G 7L + we have | YHL- n ( x i ~ (Tx)i)\ < 2, since during one iteration of 
the map at most one particle can enter the interval of sites from — n to m (from behind) and at 
most one particle can leave this interval. 

By the definition of the lower density there is a sequence of pairs (nj,rrij) 3 — °* (oo, oo) such 
that 

Xi — ► p-{x, 1). 



n, + mo + 1 . 

3 l =- n 3 



On the other hand, since | J2iJ~ n j x i ~ J2iJ- nj X^ x )i\ — ^ ; we deduce that p-(x, 1) is a limit 
point for partial sums for the sequence Tx. Therefore we need to show only that this is indeed 
the lower limit. Assume, on the contrary, that there is another limit point, call it £, for the 
partial sums for Tx such that £ < P-(x, 1). Doing the same operations with the partial sums 
for Tx converging to £ we can show that this value is also a limit point for the partial sums for 
the sequence x, and, hence, £ cannot be smaller than p~(x, 1). 
The proof for the upper density follows from the similar argument. 

By Lemma ^1] we have p±(x, 1) = 1 — p^(x, 0), which proves the preservation of the density 
of zeros as well. 

To prove that all other statistics are not preserved under dynamics we need to study it in more 
detail. Therefore we postpone the continuation of the proof till the end of this Section. 



Lemma 2.3 T* = T_ 



Proof. The action of the map T on binary configuration is equivalent to the exchange of any 
pair 10 to 01. Since the dual map describes the dynamics of holes it corresponds in this case to 
the exchange of pairs 01 to 10, which proves the statement. □ 



By a cluster (of particles) in a binary configuration idlwe mean a collection of consecutive 
positions x m ,x m+ i, . . . ,x n such that n — m > 1, X{ = Vi 6 (m, . . . , n) and x m -\ = x n+ \ = 0. 
After each iteration of the map T the last particle in the cluster moves away (i.e., (Tx) n = 0) and 
either appears a new element in the cluster from the left (i.e., (Ta;) m _i = 1 and (Tx) m -2 = 0), 
or the first particle preserves its position m. Therefore the number of particles in a given cluster 
cannot increase, and the time up to the moment when the cluster length shrinks to 1 (i.e. it 
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disappears) we shall call the life-time of the cluster or the number of iterations which are needed 
for a given cluster to disappear. 
Define an integer-valued function 

I(x, i) := max{fc < i : p(x[k, i], 1) = p(x[k, i], 0)}. (2-1) 



Lemma 2.4 Let Y17=m x i = n — m + 1 (i.e. the positions from m to n correspond to a cluster 
of particles and let I(x,n) > — oo, then after exactly ^(n — I(x,n) — 1) iterations (which is equal 
to the number of ones minus one in the word x[n — I(x,n),n]) this cluster will disappear. If 
p+(x) < 1/2 then Vi G TL we have I(x,i) > — oc. 

Proof. Let n 2n := {A G {0, l} 2n : A 2n = 1, I(A,2n) = 1}, where I(x,i) is defined by the 
relation (2.1) and thus p(A, 1) = |^4|/2. Observe that if A G 2n than for any < m < n 



and any word B G {0, l} 2m such that Bi = Ai V0 < i < 2m we have B g" $7 2m . Therefore we 
shall call the words from Q 2n minimal words (or minimal intervals) corresponding to clusters of 
particles in their ends. 
Consider a map T : J7 2n — > 7Z? n ~ 2 defined by the relation: 

(TA) i+1 = Ai + min{Ai_i, 1 - A^ - min{yli, 1 - A l+1 }. 

Observe that this is a shift to the right of the action of our map T. We shall prove that for each 
n we have r : Q 2n -> 2n " 2 . 



A= 0001011011 00100111 
T 1 A= 00011011 001011 
r 2 ,4 = 1 1 1 1 1 

r 3 ,4= ooii oi 

r 4 ^ = oi 

Examples of the action of T* on Q 10 (left) and O 8 (right). 



Let A 6 Vt 2n and let £4 be the position of the last in A. For each word A G £l 2n define a new 
word A 1 G 2Z 2n ~ 2 as follows: 



A' 



Ai if i < 2n — 2 and i ^ Qa 
1 if i = (a- 



Then \A'\ = \A\ - 2 = 2n - 2, p{A' , 1) = 1/2 and thus A' G Q 2n ~ 2 since otherwise A would be 
not minimal as well. Now using the following simple identity: 




if a -^{1,2} 

if i = Ca - 2 

if i = a - 1 

if i = \A\ - 2, 



we get TA G VL 2n ~ 2 . Example: A = 001011, A' = 0011. 

It remains to show that if we have a cluster of particles located in the end of a minimal 
configuration A G VL 2n then this cluster (i.e. particles at sites from Qa + 1 to 2n) will vanish after 
n — 1 iterations. Observe that after one iteration of the map T the cluster either preserves its 
length, or the length decreases by one (when two positions immediately preceeding the cluster 
are occupied by two zeros). The map T defined above controls this process since for each A G f2 2n 
and each < t < n — 1 the last positions starting from (Cr*A + 1) correspond to the cluster 
under study. □ 
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Corollary 2.5 The dynamics of a cluster of particles depends only on the distribution of parti- 
cles lying below the cluster. Moreover, for a given cluster of particles only particles belonging to 
its minimal word can join the cluster. Thus, if a particle does not belong to any minimal word, 
then for each t G 5Z + U {0} its local velocity in the configuration T l x is equal to 1, i.e. it moves 
freely. 



Lemma 2.6 Let A := x[n,m] and A' := x[n',m'] be two minimal words in the configuration 
x £ X. Then the inequality m < ml yields either m < n' (i.e. A n A 1 = $), or n' < n (i.e. 
A c A'). 

Proof. Assume on the contrary that n < n' < m. Then by the definition of a minimal word 
we have 

— {m! — n' + 1) = p(x[n' , m'], 1) • {m — n' + 1) 

= p(x[n' , m], 1) • (m — n' + 1) + p(x[m, m ], 1) • (rri — m + 1) 
> — (m — n + 1 + m — m + 1) = — (m — n + 2). 

We came to a contradiction. □ 



Corollary 2.7 For a given particle £ in a configuration x G X let the length of the largest 
minimal word to which £ belongs is 2n. Then for any t > n — 1 the local velocity of the particle 
£ in the configuration T l x is equal to 1. 

Observe that in the case v = M = 1 the set Free = {x G X\ : XiXi + \ = Vi} and is 
the union of 'free' particles (i.e. particles having velocity 1), while its dual Free* = {x G 
X\ : (1 — Xi){\ — Xi+i) = Vz} corresponds to 'free' holes (i.e. to holes having velocity —1). 

Lemma 2.8 Let p+(x, 1) < 1/2, then each cluster of particles in the configuration x will disap- 
pear after a finite number of iterations and dist(T*x, Free) *z^5 > o. If there is a cluster of particles 
having an infinite minimal word (i.e. which does not vanish in finite time), then p+{x, 1) > 1/2. 
If p+{x, 1) > 1/2 then there are clusters of particles with arbitrary large (but may be finite) 
minimal words. 

Proof. Let x G X satisfies the assumption that p + (x,l) < 1/2 and let the segment x[n,m] 
be a cluster of particles. Then there exist a pair of integers n', m' such that infty < n' < n < 
m < m' < oo and x[n',m'] is the largest minimal word covering the cluster of particles x[n,m] 



(otherwise this would contradict to the definition of the upper density). Hence by Lemma 2.4 
after at most (m' — n') /2 iterations this cluster will disappear and all particles will become free. 
Since this argument can be applied to any cluster of particles, this yields the first statement. 

Assume now that the minimal word of a cluster of particles x[n,m] is not bounded. Then for 
any n' < n we have x[n',m] > 1/2 and thus for any k G 7L + we have 

p(x\n — k , m + k , 11 > -5 p i n — k , m , 1) 

m — n + + 1 

m-n + k 2 + 11 / k \ 1 k^oo 1 

> : r- =1 



m + k-n + k 2 + 12 \ m + k-n + k 2 + 1)2 2 

Therefore p+(x, 1) > 1/2, which proves the second statement. 
The last statement is an immediate consequence of the definition of the minimal word. □ 
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Lemma 2.9 Let p_(x, 1) > 1/2 then dist(T**x*, Free) ^ 0. 



Proof. By Lemma ^J] we have p+(x*, 1) = 1 — p~{x, 1) < 1/2. On the other hand, Lemma 2.5 
shows that asymptotic properties of the maps T and T* coincide, thus we can apply the state- 
ments of Lemma 2^8 for the case of (T*)*x* to prove the desired result. □ 



Note now that there are configurations not satisfying the assumptions of Lemmata 2.9,2.8 
which still converge to Free U Free* under the action of the map T. Indeed, let y = . . . 111000 . . . 
and let the index correspond to the first in y. Observe that p~(y, 1) = < p+(y, 1) = 1, 
however dist(T*y, Free) = 2~(* +1 )(l + 2~ 2 + 2" 4 + . . .) = § • 2~* 0, since for large t the 
'central' part of T l y will be occupied only by free particles. On the other hand, for y* = 
. ..000111. .. we have dist(r*y*, Free) = T x + 2" 3 + . . . + 2" 2n+1 + ... = 1/3 for each t G E+ 
while dist(T*V, Free) = § • 2~* *-=^ 0. 

Lemma 2.10 For any x G X we have dist(T*x, Free) i -^5 > if and only if p + (x[l, oo], 1) < 1/2. 
If additionally there exist a pair n,m G 7L + such that p(x[i,i + m — 1], 1) < 1/2 for each i > n 
we have dist(T*x, Free) < Const 2 . 

Proof. Observe that p+(x[l, oo], 1) < 1/2 implies that there exists N G 2Z such that the 
life-time for each cluster of particles lying to the right from iV is finite. On the other hand, the 
distance from the most right cluster of particles in T l x located to the left of iV to the position 
N grows with t linearly. This proves the first statement and shows that the rate of convergence 
might be smaller than 2~* only if the life-time of clusters of particles lying to the right from a 
sufficiently large position N is not bounded. The additional assumption guaranteers that this 
cannot happen, which yields the second statement. □ 

Let M(X) be the set of probabilistic translation invariant measures on X and let /i[^(x)] := 
/ (p(x) dp(x) for p G M(X), in particular, p[xo] := p(x G X : xq = 1). Consider a subset 
of Ai(X) corresponding to measures in the space of sequences with weak dependence between 
coordinates (exponentially fast decay of correlations): 

M P {X) := {p G M{X) : p[x ] = p, \p[x ■ x k ] - p 2 [x }\ < Ce^}. 

for some C,a and VA; G 2Z. Note that, e.g., a product (Bernoulli) measure p p G M P (X). 

Lemma 2.11 For any p G Ai p (X) we have p(x, 1) = p for p-a.a. x G X , and thus p(x G 
X : p_(x) < 1/2 < p + (x)) = 0. 

Proof. Let S n>m (x) := J2iL- n x i- Then p[S n>m (x)] = p and by Chebyshev inequality Ve > 
we have 

p{x G X : \S n . m (x) -p\ > e) < • p[(S n . m (x) - p) 2 }. 

e z 

A straightforward calculation shows that p[(S n ^ m {x) — p) 2 ] < n _Sl +l and thus 



p x G X : 



1 "' 



n + m + 



Xi-p 



n+m+1 

> e \ < — ;r — > 0. 

[n + m + lje 2 



fit- ^oo 

Therefore p±(x[—n,m],l) ' — > p in probability, which yields the existence of the density 
p(x, 1) = p for p-a.a. x G X and thus the statement under question. □ 

Denote by Per n (T) := {x G X : T n x = x} the set of n-periodic (in time) trajectories of the 
map T and by B(Y) := U n >oT~ n Y the basin of attraction of a subset Y C X. 
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Lemma 2.12 T : X — > X is a Lipschitz continuous map in the topology induced by the metrics 
dist(-, •). For each n G 7L + there exists an n-periodic trajectory, and all periodic trajectories are 
unstable. Clos(£(Free U Free*)) = X, (Free U Free*) n Peri(T) = 0, and fi p (B (Per ^T))) = 
while Clos(B(Peri(T))) = X. 

Proof. Let us start with the Lipschitz continuity. Consider two configurations i/i/£l and 
assume that — n < is the largest negative index and m > is the smallest nonnegative index 
of sites, where they differ, i.e. for all —n < i < m we have Xi = y%. Then we have 

2 _ n + 2 _ m < digt ^^ ^ < 2 ^_„, + 2 _ m ^ 

On the other hand, a straightforward calculation shows that the closest to the origin from the left 
side differing coordinates of the configurations Tx and Ty belong to the set {— (n+1), —n, — (n — 
1)}, while the closest from the right side belong to {m — l,m,m+ 1}. Thus 

2 -(n+i) + 2 -(m+i) < dist(Tz,Ty) < 2 (2-("- 1 ) +2"( m - 1 )). 

Therefore 

I ^ 2 ~( n+1 ) +2-( m+1 ) dist(Tx,Ty) < 2(2~( n - 1 ) + 2-( m ~ 1 )) _ 
4~ 2(2-™ + 2~ m ) ~ dist(x,y) ~ 2~ n + 2~ m ~ 

For a given n G Zj + consider a space-periodic configuration x £ X with the (space) period n, 
e.g. Xj = Xi +n Mi. Then it is immediate to show that for any t &7L the configuration Tx is again 
space periodic with the same period n and converges either to Free, or to Free*, depending on 
its density. This gives a construction of the n-periodic (in time) trajectories. 

The structure on the set of fixed points Peri(T) is a bit more involved: 



Peri(T) := \ x (n) G X : 



(n) _ J if i < n 
1 I 1 otherwise 



Indeed, assume that Tx = x, then either x does not have zero coordinates, or all coordinates 
starting from, say, n-th, should be equal to one. Now for G Peri(T) we define y( n ' m ) G X 
such that 

(n,m) _ J if i < n, or i > m 
1 1 otherwise 



Vi 



for some m > n. Then dist(x( n ), y( n ' m )) = 2~ m 0, while dist(TM n ), T t y( n > m '>) *-=^ 

2-(«-i) _^ q xhus for each e > there is configuration x' = x'(e) such that dist^™), x') < e 
and T l x' -f* x^ as t — » oo, which yields instability. 

Observe now that the set K := X \ (#(Free U Free*)) = {i£l: p_(x) < 1/2 < p+(x)} has 
/ip-measure zero, since for each ^,-typical trajectory the lower and upper densities coincide. 

Consider now an arbitrary configuration x € X and a sequence of configurations {y^ n '} n defined 

as = |^ ^enrise • Then dist(x,y(")) < 2~ n+1 , on the other hand, y^ Peri(T), 
which proves the last statement. □ 



Remark. In the case of a finite cluster of particles its last particle immediately leaves under 
the dynamics. This is not the case for clusters not bounded from the right, which explains the 
existence of fixed points. 

We shall say that a closed T-invariant set Y is a weak attractor if /j, tc {(B(Y)) > 0. A weak 
attractor Y is called a Milnor attractor if /x re f (B(Y)\B(Y')) > for any proper compact invariant 
subset Y' C Y (see, e.g. @). 



9 



Lemma 2.13 The set Free U Free* is an week attractor with respect p re f = p p , but not a Milnor 
one, moreover it is not a topological attractor. 



Proof. The sets Free and Free* are closed, since they contain all their limit points. Let 
ZW ■= {x £ Free : p(x) = p} - this is an invariant set and p p (B(Z^)) = 1. Denote now 
by a single configuration from together with all its left and right shifts. Clearly 

fi p (Z(P) \ (Z( pS) \ Z'( p ))) = u p (Z'( p ^) = 0. Observe that the points from the complement to the 
basins of attraction of Free and Free* are everywhere dense, which proves the absense of included 
open sets. The last statement follows from the fact that the basin of attraction does not contain 
any open set. □ 



Proof of Theorem 1.1 (continuation). Let us prove now that for any word A with \A\ > 1 the 
density p±(x, A) is not preserved under dynamics. There might be 3 possibilities: p(A, 1) < 1/2, 
p(A,l) > 1/2 and p(A,l) = 1/2. We start from the first case. Clearly p(A,l) < 1/2 yields 
p(A, 00) > 0. Consider a configuration x := (A 1 1 . . . 1 ) , where x = (B) = . . . BBB . . . stays 

2\A\ 

for a space-periodic configuration. By the construction p(x, 1) = (p(A, 1) • \A\ + 2|j4|)/(3|A|) > 
2/3 > 1/2. Therefore T l x Free* and hence p(T t x, 00) *z^5 > q Assume now that the density 



is preserved, i.e. p(T t x,A) = p(x,A) Vi. Then by Lemma |2,2| we have 

p{T*x, 00) > p(T f x, A) ■ p(A, 00) = p(x, A) ■ p(A, 00) > 0, 

while the left hand side vanishes when t — > oo. We came to the contradiction. 

If p(A, 1) > 1/2 we shall follow a similar argument, considering another space-periodic config- 
uration x := (A00...0). 

2\A\ 

In a more delicate case p(A, 1) = 1/2 we do the folowing. If additionally p(A, 11) > we follow 
the same argument as in the case p(A, 1) < 1/2 to show that p(x, 11) > 0, while p(T l x, 11) * z ^5 > 0. 

If p(A, 00) > we follow the case p(A, 1) > 1/2 to show that p(x, 00) > 0, while p(T x , 00) ^ 
0. It remains to consider the case when p(x,ll) = p(x, 00) = 0, i.e. A = 1010... 10 or 
A = 0101 . . . 01. In the first of these case we choose x := (1^40). Then 

p((lA0),A) = lim -— ? = — i — < - = lim 7 ^^ = p(T((lA0)),A). 

In the second case we choose x := (0^41) to come to a similar contradiction. □ 



Proof of Theorem 1.2 in the case v = M = 1 follows from Lemmata 2.8-2.12. □ 



Proof of Theorem L3 in the case v = M = 1. We have the following identity: p±(x, 1) = 
p±(x,10) +p±(x, 11). Hp + (x,l) < 1/2 thenT*a; -> Free and p + (T*x,ll) -> 0, thus $±(T t x) = 
p±(T t x, 10) = p±(T t x, 1) - p±{T t x, 11) -» p±{T t x, 1) = p±(x, 1). The situation 1) > 1/2 

can be reduced to the previous one by going to the dual configuration. 

Consider now the case p~(x, 1) < 1/2 < p+(x, 1). By definition there exists a sequence of pairs 
of positive integers ra-,m- — > oo such that p(x[—n' i ,m' i \,l) * z ^> p_(x,l) < 1/2. For each i we 
choose integers ni > n^, iri{ > m! i to be the smallest integers satisfying the condition that — — 1 
is the ending point and + 1 is the starting point of some nonoverlapping minimal intervals 
of the configuration x. If there are no more nonoverlapping minimal intervals in the considered 
direction or the segment x[— n'^m^] intersects with an infinitely long minimal interval we set 
ni := n\ or rrtj := m\ respectively, depending on the direction where this event occurs. Clearly, 
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we have p(x[— ni,m,i], 1) < p(x[—n' i ,m / i ],l) and thus p(x[— rii, mil, 1) * :: ^5 > p_(ar, 1) < 1/2. By 
the definition of minimal intervals after tj := (n« + mj)/2 + 1 iterations all clusters of particles 
inside of the segment x[— n^m^] will disappear and all particles will become free. Therefore we 
can again apply the same argument as in the case p+(x, 1) < 1/2 and obtain the relation for the 
lower limit of the flux. To obtain the relation for the upper limit one should consider the dual 
configuration. □ 



3 Dynamics of fast particles (T v , X) 

Note that the analysis of dynamics of the slow particles model (T, X) is divided logically into 
two parts: first, we study low density initial configurations x € X with p+(x, 1) < 1/2, and 
then for high density configurations x £ X with p~(x, 1) > 1/2 we pass to the dual ones using 
the property that p+{x*, 1) < 1 — p~(x, 1) < 1/2 and argue that the dual map T* = T_i has 
exactly the same asymptotic properties as T. The problem with the fast particles model (T v , X) 
is that the dual map map T* ^ T- v in this case, and, in fact, has a very nontrivial dynamics. 
Namely, T* v corresponds to the situation, known in physical literature (in the case v = 2) as 
a traffic model with 'smart drivers', who anticipating the motion of at most v cars ahead, may 
move to an occupied site ahead of it with the maximal velocity 1. Example for the case v = 2: 
(OHIO) (01011). 

Therefore since we are unable to study directly the dual map in this case and according to the 
entire ideology of this paper, we elaborated a reduction to the main case v = M = 1 based on 
the following consideration. Note that under the action of the map T on x G X each pair 10 
goes to 01 (i.e. the position of a particle and a hole are exchanged). Therefore T is equivalent 
to the substitution rule 10 — > 01. To apply this idea to the case of T v we introduce an alphabet 
A v := {0i, O2, • • • , V , 1} with v + 1 symbols and a map C v : X — > = X v defined as follows: 
for each segment x[i,i + n + 1] = 10^01, we set C v x[i,i + n + 1] := 10„ . . .0„ ^ n -\_ n /v\v^-- If 

\n/v\ 

n — [n/v\ v = we shall drop the last element in C v x[i,i + n + 1]. It remains to define the 
action of C v on 'tails' of x consisting of only zeros, which we set according to the following rules: 

. . . 0001 ... ... 0„0„0„1 . . . and . . . 1000 . . . . . . lO^C . . .. 

Now we are ready to define the substitution map S v : X v — > X v acting in the set X v according 
to the set of v substitution rules 10j — > Ojl for < i < v, which generalizes the substitution rule 
for the slow particles dynamics for the case of v different types of holes. 

To study the life-time of clusters of particles in configurations x £ X v we introduce also a new 
map T := C V C V ^S^Cv. 

Lemma 3.1 T v = C~ 1 S V C V , and T™ = C~ l f n for any n G ZZ+. 

Proof. Straightforward. Note only that the map C V C~ X needs not to be identical, example: 

CC-^lO^-iOil) = (b%l). 

□ 

Observe that the map T acts on X v in exactly the same way as T acts on the space of binary 
sequences, namely T moves each particle by one position forward if there is no particle there 
or the particle preserves its position otherwise. So the only difference is that now we have v 
different types of zeros, instead of the only one type in the case v = 1. 
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Therefore to study the life-time of a cluster of particles we apply a similar machinery as in the 
case of the slow particles. Denote 



Ind(a) ■= \ . V _ i > : = max{£; < % ; ^ Ind(Aj) < 0}, 

l z 1±a - u » j=k+1 

Q n := {A £ A n M : A n = l, I(A,n) = l}, 

where Aj is (as usual) the j-th element of the word A. Consider a map T defined on words of 
length n 6 7L + as follows: (TA)i := {TA)i + \ for all i = 1, 2, . . . , n — 1. 

Lemma 3.2 T : S7 n — > il n_2_ ^ ; where < £ < n — 1. JTie life-time of the cluster of particles in 
the end of a word A E fi n is e^ua/ to (p(j4, 1) • \ A\ — 1). 



Proof follows from the same argument as the one of Lemma 2.4. The only difference is that 
due to the action of V the number of elements in TA may become smaller than \ A\ — 2, since the 
action of C V C~ X may decrease the number of 0j. On the other hand, during the one iteration of 
the map T only one elemement 1 disappears from (the right hand side) of A, i.e. p(A, 1) • \A\ = 
p(TA, 1) • \TA\ + 1. Therefore the number of iterations needed for the cluster of particles in the 
end of the word A to disappear is equal to the number of ones in the word A minus one. □ 



A= 2 2 0i 1 0i 1 0i 1 1 

F l A= 2 2 1 0i 1 1 

T 2 A= 0i 2 1 1 

F 3 A = 2 1 

Example of the action of T* on ft 9 with v = 2. 



Proof of Theorems [lTH]l.2|JL~3l for the case v > 1, M = 1 follows immediately from Lem- 



mata 3.1,3.2 and the reduction to the case v = 1 obtained there. □ 



Consider now a special case of superfast particles corresponding to the choice of maximal 
velocity v = oo. Denote 

X (oo) :={xeX : Vneffi 3m, m > \n\ : x m = 1, x_ m ^ = 0}, 

i.e. the set of binary configurations having no infinitely long right 'tails' of zeros or left 'tails' 
of ones. Then the maps T^^T^ : — > are well defined. The substitution rule 

1(L^C)1 — > 10jl Vi € 7L + maps X^°°^ — > Xoo. Strictly speaking, the latter has an infinite 

i 

alphabet, however all arguments applied in the case of finite v work as well. Moreover here the 
situation is even simpler, because between each pair of consecutive ones there is only one zero 
with a certain finite index: . . . lOjlOj .... Thus the dynamics of (T^, Xqo) is equivalent to the 
dynamics of free particles, which gives the flux &±(x) = 1 — p±(x, 1). 

4 Dynamics of multi lane flows (T Vj m, Xm)- 

Reduction of Ti m to the direct product of M maps T 

The model of a muti lane flow of slow particles on a finite lattice has been introduced in [|TJ] 
and generalized for the case of an infinite lattice TL in M , where statistical properties of regular 
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Figure 2: 'Sawtooth redirection': (a) action of the map Si on individual particles marked by 
squares with their relative numbers (M = 3); (b) symbolical representation of Six, - one 

of intervals of monotonicity. 



initial configurations have been obtained. However the approach used in |H] does not allow to 
study the dynamics of general initial configurations, which we shall consider in this Section using 
a completely different method. 

Our first aim is to redistribute a configuration x G Xm into M binary configurations {x^ G 
X = Xx}f =l , such that T\ M x = J2j T txij) for a11 t£7L+\J {0}, where the notation x = ^ 
means that xi = J2j x i f° r each i G 7L. To solve this problem we introduce a sawtooth 
redirection Si : Xm —> {X\) M with Six = {x^'}jL 1 of a configuration x € Xm to a collection of 
binary configurations {x^}jL 1 with the starting point at site I G 2Z: 




1 if % > I and j G (0^ 0| = ; a? fc ] 

1 if i < I and j G (0^(-x fc ), 0^ (-**)] 

otherwise, 



where a © 6 := (a + 6 — l)(modM) + 1 and 0^L n xi := x„ © . . . © x m . In other words, for the 
configuration x G Xj^ we construct a beinfinite 'staircase' starting from the site / with the i-ih 
stair of height Xi and then redistribute the result modulo M (preserving the site number) among 
M binary configurations {x^'}j^ =1 . 

With some abuse of notation we shall refer to (Six)^ = x^ as the j-th lane of Six = {x^}jL 1 
and denote the action of the direct product of maps T v applied at Six as T v Six := {T v x^}jL 1 . 

Example for the case v = 1, M = 3 and the starting site / defined by the relation x [I — 1 , 1 + 1] = 
211: 

...10100110... ...01010101... 

Si(... 11211221...) = ...00101010... ...* 0010101... = S l+1 (. ..* 112121 *.. .), 

...01010101... ...* 010101*... 

where the unknown positions are marked by *. 

Symbolically the sawtooth redirection is shown in Fig. ||(b) by curvilinear lines corresponding 
to sawtooth rows of ones, open circles mark the intersections of these lines with the 'lanes' j, j', 
i.e. the positions where x^> or x^ are equal to 1 (all other positions on these lanes are occupied 
by zeros). 

Theorem 4.1 For any x G Xm and I G TL and Si(x) = {x^^'Li we have 

(a) x = Z j (S l x)W, 

(b) \p(xW)[n + l,n + fc],l) - p(x^[n + l,n + k], 1)| < 1/k Vj,j' G {1,...,M}, n G 7L and 
k G 7L + , 

(c) Si +k x = {(Six)^}f =1 Vk G 7L + , 



13 



(d) Vi> > 1 we have T v Six = Si+t-x' for some £ € {0,1,... ,v} and x' € Xm which doesn't 
depend on I, 

(e) TSix = Si + ^(Ti >M x) for some £ € {0, 1}. 

Proof. The statement (a) follows immediately from the definition of the sawtooth redirection, 
because during the redirection each particle preserves its position i. 
The property (b) is equivalent to the assumption that 



i=l i=l 



< 1, 



i.e. that the number of particles in the same segment of different 'lanes' j,f can differ at 
most by 1. According to the 'sawtooth redirection' for any given finite segment of integers 
re + 1, n + 2, . . . ,n + k the number of intersections of the curvilinear lines in Fig. §(b) with the 
horizontal line at height j differs from number of intersections with the horizontal line at height 
j at most by one. This immediately yields the property (b). 

The collection of binary configurations S[(x) has a row of ones at site i of height k if and only 
if Xi = k, and the change of the starting point I of the redirection only changes cyclically the 
starting point 1 of the enumeration of lanes x^\ This proves the property (c). 

Observe now that the definition of Si(x) is equivalent to the existence of a partition of TL into 
segments such that x\_ = 1, = 1 (except for the most left segment where i_ = — oo 

and the most right one where i+ = oo) and for any 1 < j < M there exists the only one 
i G [£-,*+] such that x*f' = 1. Indeed, according to the 'sawtooth redirection' the curvilinear 
lines in Fig. ||(b) have the property that the intersection with the horizontal line at height j 
occurs not earlier than with the horizontal line at height j' > j (the curvilinear lines may have 
vertical segments). To simplify the notation we shall say that Si(x) is monotonous on 

Consider the interval of monotonicuty which starts from I, i.e. i_ = I. We set £ to be 

equal to the minimum of the number of not occupied positions in x^ ahead of the site i- (which 
is occupied by 1). Then under the action of T v the particle at the site i_ of the 1-st lane moves 

by £ positions to the right. Observe that all particles on the other lanes in the segment [z , i-)-] 

have at least £ not occupied positions ahead of them, and therefore all these particles will move 
at least £ positions to the right. Thus to prove that the monotonicity is preserved it is enough 
to note that the particle on the lane M cannot move further to the right than the first particle 
on the first lane of the next interval of monotonicity. Indeed, the latter is a trivial consequence 
of the definition of intervals of monotonicity. This finishes the proof of the statement (d) except 
the last part, which follows from the statement (c). 

To prove the statement (e), observe that by the definition of the map T\^m (see Section ||) a 
particle at the site i of the lane j can switch to the lane j' if and only if x^[i,i + 1] = 11 and 
x^ + = 00, which contradicts to the definition of the intervals of monotonicity. Therefore 
under the sawtooth redirection no particle in Si(x) will change its lane. □ 



Corollary 4.1 The sawtooth redirection gives a simple constructive way to rearrange vehicles 
in a multi lane traffic flow between lanes (preserving their positions in the flow) to achieve the 
maximal available flux. 

According to Theorem [4.1| ,(d) the map x — > J2j T v (Six)^ is well defined as a map from Xm 
into itself and does not depend on the choice of the starting site / S TL. Moreover, it can be 



shown that this formula coincides with (1.1) in the case v = 1, and it clearly coincides with T v 
in the case M = 1. Therefore we use this relation as a definition of the dynamics of a general 
multi lane flow in the case v, M > 1, namely we set T Vj m% '■= J2j T v (Sqx)^' . 
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Proof of Theorem 1.1 for the case v,M > 1 and A C X follows now from the sawtooth 
redirection, which gives the reduction to the one-lane case. It remains to show that the statistics 
of more general words A C Xm with \A\ = 1 might be not preserved under dynamics. The reason 
of this is that if M > 1 the multiplicities might be not preserved. Indeed, let a E Am \ {0, 1}. 
Then p((a(M - a + 1)0), a) = ±(1 + |_(Af - a + l)/aj +0), while 

p{T 1M (a(M - a + 1)0), a) = p«l(a - 1)(M - a + 1)), a) 

= \(0 + + L(M - a + l)/aj + 0) < p((a(M - a + 1)0), a). 

□ 



Proof of Theorems |L2]]L3| for the case u, M > 1. Consider a configuration x E X^f . According 
to Theorem [OJ,(b) for Sqx = {x^}f£-y we have Vn,m E K + that 

|p(a?^[— n, ml, 1) — p(x^ H— n,m], 1)1 < . 

m + n + 1 

Thus going to the limit as n, m — > oo and using Theorem |4.l|, (a) we get p±(x^\ 1) = jjp±(x, 1) 
for each j E {1, . . . , M}. Therefore the application of the results obtained in Sections |||| in the 
case of one-lane flows (i.e. in the case of the map T v ) proves the statements under question. □ 



5 Rate of convergence: (space) periodic, regular, and typical 
initial configurations 

In this section we study the rate of convergence of various statistics to the corresponding limit 



values whose existence have been established in Theorems 1.3,1.2, Since we have shown that 



the analysis of T v m in all cases can be reduced to the case of v = M = 1, we consider in this 



section only the dynamics of slow particles and shall consider the proof of Theorem 1.4 only for 
this case. 

We start with periodic in space configurations. Clearly each n-periodic in space configuration 
x E X can be represented in the form x := {A) with a binary word A of length \A\ = n. 

Lemma 5.1 Let x := (A) with \A\ = n, then p±(x,a) = p(A,a) for a E {0,1}. The space of 
n-periodic in space configurations is invariant under the action of the map T and after at most 
|_n/2j +1 iterations any configuration from this space belongs to Free U Free*. 

Proof. Straightforward. □ 

The only nontrivial question related to this space is the length of the transient period for a 
given configuration 

Lemma 5.2 Let x := {A) with \A\ = n and p(A, 1) < 1/2 and let B be the longest minimal 
word in A. Then the length of the transient perion is equal to \B\J2 — 1. 



Proof. This is an immediate consequence of Lemma |2.4| . □ 



Consider now a generalization of the space of periodic in space configurations - the space of 
regular configurations, proposed in [Q, This space is defined as follows: 

Reg(r, ip) := {x E X : \p(x[— n, m], 1) — r| < ip(n + m) Vn,m£ffi + }, 
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where r G [0, 1] is a constant, and the nonnegative function ip(n) — > is assumed to be a 
strictly decreasing. 

Lemma 5.3 Let x G Reg(r, ip), then the density p(x, 1) := r is well defined, and if ^ 1/2 then 
there is a constant r such that the life-time of any cluster of particles in x does not exceed r. 

Proof. First, observe that if x G Reg(r, -0) then we have 

^^ 1 ) = n ^ r 5oo( in? Hl) = r, 

and thus p(x, 1) = r is well defined. Assume now that x G Reg(r, ip) with r < 1/2, then, since 
we deduce that there is a positive integer r such that <P(t) < 1/2 — r. Then we have 
p(x[n, n + t], 1) < 1/2 for any n € 7L, which yields the claim of the lemma due to Lemma |~~ 



The case r > 1/2 follows from the same argument but applying Lemma 2.£ instead. □ 

Now we proceed to study more general initial configurations. 

Lemma 5.4 Let x G X satisfies the assumption that there exists a number 7 G (0, 1) such 
that Vn G 2Z+ and any word A C x[—n,n] with \A\ > 2jn we have p(A, 1) < 1/2. T/ien 
dist(T*x, Free) < 2~*/ 7+1 /or any t G 2Z + . 7/ x* satisfies the same assumption, then we have 
dist(T*x,Free*) < 2- l ^ +l . 

Proof. Consider only those n E 7L + for which the largest minimal words containing a cluster 
of particles in x{— n,n] also belong to x[— n,n]. By the assumption of Lemma the length of 
the largest minimal interval containing in the segment x[—n,n] does not exceed 2777,. Therefore 
the corresponding clusters of particles with disappear after at most 771 iterations, and thus 
for all sufficiently large t € 7L + all particles in the segment T t x[—t/j,t/'y] will become free. 
Thus the closest to the origin nonfree particle can appear not earlier as at site i/7, which gives 
the desired estimate of the rate of convergence. The second statement follows from the same 
argument applied to the dual map. □ 



Lemma 5.5 Let x £ X satisfies the same assumption as in Lemma 5.4, then 

1 n 

hmsup— V(T n x,l) = F 1A (p(x,l)). 



Proof. Observe that ^ E?=- n v ( x > x ) = p( x 

[—n,n], 10). Applying the same argument as in 
the proof of Lemma we see that after n iterations the segment T n x[— n/7, n/7] contains only 
free particles. Therefore p(x[—n,n], 10) = p(x[—n,n], 1), which yields the desired equality. □ 



Corollary 5.6 The statements of Lemmata 5.4,5.1 remain valid if instead of Vn € 7L + we 



assume that n belongs to the subset of 7L + of density 1. 



Lemma 5.7 V7 € (0,1) for p p -a.a. configurations x G X the set of n G for which any 

word A C x[—n,n] with \A\ > 27n satisfies the inequality p(A, 1) < 1/2, has the density 1. 
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Proof. [] Let {xi}??^ be a Bernoulli sequence with the density V(x{ = 1) = p < 1/2 for 

all i G 2Z. Introduce a sequence of functions y n {r) := 2n +i Y^i^-n^^ x * depending on a real 
variable r G [7, 1], and consider a functional 

4>{y( T )) : = su p su p -(y( r + s ) - v( T )) 

re[0,l-7] 7<s<1-t S 

defined in Skorohod space of functions y(r). Then the quantity under question is the probability 
V{4>{y n {j)) < 1/2 Vr G [7,1])- Since y n {r) converges in probability for a given r to y(r) : = 
pr and the functional (f) is continuous, 0(y n (r)) converges to (j)(y(r)) (functional law of large 
numbers). Thus we have 

V{4>{y n {r)) < 1/2 Vr G [7, 1]) -> 7>(0(y(l)) < 1/2) = 1, 

where the rate of convergence 

(P(<KVnQ)) > l/2)) 1/n "-=3? VM 1 ~P) follows by the combina- 
tion of the large deviation principle for the functions y n ( T ) an d the contraction principle (see, 
e.g., 0). □ 



Corollary 5.8 Results of Lemmata 5.4.5.5.5.7 prove Theorem \1.4\ in the case v = M = 1. 



6 Dynamics of measures and chaoticity 

In this section we shall study the action of the map T v m m the space A4(Xm) of probabilistic 
measures on Xm- This action is defined as follows: T Vi m[J<(Y) := ^(T^^Y) for a measure 
H G M(Xm) and a measurable subset Y C Xm- A measure \x G A^(Xm) is called translation 
invariant if it is invariant with respect to the action of the shift map a : Xm — * Xm- 

Lemma 6.1 If \i G M{Xm) is translation invariant then this property holds for M \x Vf G 

Proof. We have T^^Y) = (j,(T~ M Y) = ^aT^Y) = ^(T~ M aY) = T^^aY). □ 

One might expect that under the action of the map T v m any translation invariant measure 
should converge to a Bernoulli one. Indeed, 

7> p (x G X : x = l) = fip{x G X : x[0, 1} = 11) + fi p {x G X : x[-l, 0} = 10) 

= n P {x G X : x[0, 1] = 11) + fip(x G X : x[0, 1] = 10) = fi p (x G X : x = 1) 

On the other hand, the product structure is not preserved even in the case of the model of slow 
particles. 

Lemma 6.2 The measure Tfj, p is not a product one for any < p < 1. 
Proof. We have 

7> p (x G X : x[0, 1] = 11) = fj, p (x G X : x[0, 2] = 111) + fj, p (x G X : x[-l, 2] = 1011) 

= p 3 +p 3 (l-p) =p 3 (2-p) ^p 2 = fj, p (x G X : z[0,l] = 11). 

Thus the measure T\i v does not have the product structure. □ 



The idea of this construction, based on the large deviation principle, was proposed by A. Puhal'skii. 
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It is of interest that in the case v > 1 even the average value p v {x G X : xq = 1) is not 
preserved under dynamics. Indeed, 

V 

T v p p (x G X : x = l) = p p (x G X : x[0, 1] = 11) + ^ fi p (x G X : z[-l, 0] = 10.) 

i=l 

v 

= p 2 +p(i -p) + . . . +p(i -p) v =p + p^(i - P y 

i=2 

> P = Hp(x G X : xq = 1). 

Note that in the case of the slow particles model (Ti t \,X) some results about the set of T^i- 
invariant measures and mathematical expectations of the limit flux with respect to them were 
studied in |IJ. 

In Q it has been proven that the dynamical system (Ti t M, Xm) is chaotic in the sense that its 
topological entropy is positive. Moreover this paper gives an asymptotically exact (as M — » oo) 
representation for the entropy. The extension of this result to the case {T V) m , Xm) with v > 1 
is straightforward. 

7 Passive tracer in the 1-lane flow of fast particles 

Let T*x, v > 1 describes the 1-lane flow of particles and let at time t the passive tracer occupies 
the position i. Then before the next time step of the model of the flow the tracer moves in 
its chosen direction to the closest (in this direction) position of a particle of the configuration 
T^x. For example, if the going forward tracer occupies the position 2 and the closest particle in 
this direction occupies the position 5, then the tracer moves to the position 5. Then the next 
iteration of the flow occurs, the tracer moves to its new position, etc. 

To be precise let us fixed a configuration x G X with p-(x, 1) > and introduce the maps 
t^t :7L — > 7L defined as follows: 

T^i := min{j : i < j, Xj = 1}, T~i := max{j : i > j, Xj = 1}. 

Then the simultaneous dynamics of the configuration of particles (describing the flow) and the 
tracer is defined by the skew product of two maps - the map T v and one of the maps r. ± , i.e. 

(x,i) -+T±(x,i) := (T v x,r±i), 

acting on the extended phase space X x 7L. The sign + or — here corresponds to the motion 
along or against the flow. We define the average (in time) velocity of the tracer V(t, x) as S{t)/t, 
where S(t) denotes the total distance covered by the tracer (which starts at the site 0) up to the 
moment t with the positive sign if the tracer moves forward, and the negative sign otherwise. 

Theorem 7.1 Let x G {x G X : dist(T*x, Free U Free*) < 2 _t / 7+1 } for all t G 7L + and some 
< 7 < 1. If < p+(x, 1) < ^rj, then V(t,x) ^+ v if the tracer moves along the flow (i.e. in 

the caseT + ), andlimt^oc ^ SU ^ JV(i, x) = ^ +1 * n the opposite case. Ifp-(x,l) > 1 — ^j-j- 
and the tracer moves against the flow then V(t,x) —1. 

Remark. The assumption about the initial configurations is satisfied for p p -a.a. x G Xm (see 
Theorem |l.4|) . 

Proof. Since we assume that T*x converges to the attractor Free U Free* with the exponentially 
fast rate, then at the moment t G 7L + we have an exponentially long (in t) interval of the 
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configuration T^x consisting of only free particles or free holes (depending on the density). As 
we shall show that V(t, x) converges to a constant, then to study its value we can restrict the 
analysis to the case x € Free U Free* . 

Under the assumption < p+(x, 1) < ^j-j- we have T^x ^> Freet,. In the case of T + the 
tracer will run down one of the particles and will follow it, but cannot outstrip. Indeed after 
each iteration of the flow this free particle occurs exactly v positions aheed of the tracer. Thus 
V(t,x) *^?v. 

Consider now the case when the tracer moves backward with respect to the flow. Then each 
time when the tracer encounters a particle, on the next time step this particle moves in the 
opposite direction and does not interfere with the movement of the tracer. We assume again 
that x G Free„ and consider the case < p+(x, 1) < If on the spread of length n there are 
m particles, i.e. m obstacles for the tracer then the average velocity on this segment is equal to 
Going to the limit as n — > oo we obtain the desired estimate. 

It remains to consider the case p-(x, 1) > 1 — and thus T*x < — ^> Free*, i.e. to the flow 
where all holes move at maximal velocity —v. Thus after each iteration the tracer moves exactly 
by one position to the left (since it never can encounter a hole), which gives the limit velocity 
-1. □ 



Observe that the motion against the flow is efficient only in the case of low density of particles 
when p + (x,l) < ^j-j-. On the other hand, in the high density region in the case of the motion 
along the flow and in the region < p~(x, 1) < 1 — in the case of the motion against 
the flow the limit velocity of the tracer depends not only on the densities, but also on the fine 
structure of the configuration x. Moreover, this concerns also the case of 'untypical' initial 
configurations with < p~(x, 1) < 1/2 < p+(x, 1), when there might be arbitrary long (even 
infinite) minimal words for both particles and holes. 
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